Regular and irregular modulation of frequencies in limit cycle oscillator networks 
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Frequency modulation by perturbation is the essential trait that differentiates limit cycle oscilla- 
tors from phase oscillators. We studied networks of identical limit cycle oscillators whose frequencies 
are modulated sensitively by the change of their amplitudes, and demonstrated that the frequen- 
cies sustainably take distributed values. We observed two complex phenomena in the networks: 
stationarily distributed frequencies at a regular interval, and continuous irregular modulation of 
frequencies. In the analysis we reveal the mechanisms by which the frequencies are distributed, 
and show how the sensitive modulation of frequencies produces these complex phenomena. We also 
illustrate the topology of networks regulating the behaviors of the systems. 
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Oscillator networks are widely accepted as a standard 
model for understanding dynamical systems with many 
degree of freedom, because of the accumulation of knowl- 
edge about them [l|, 0, S 0] and the large number of 
their applications [5|. The statistical synchrony for sim- 
ple phase oscillator networks has been well investigated 
since the celebrated works by Winflee [l[ and Kuramoto 
0. Recently this field has been developed by incorpo- 
rating statistical knowledge about complex networks [6|. 
At the same time, detailed dynamic behaviors are also 
studied for oscillators showing high dimensional proper- 
ties. Various types of coherent motions are classified, and 
concepts such as phase synchronization clustering, 
and chaotic itinerancy Q have been proposed. However, 
a large number of complex phenomena still remain un- 
revealed, and the integrated understanding about them 
is still lacking compared with that about simple phase 
oscillator systems. 

Which property of a high-dimensional oscillator is es- 
sential for generating complex phenomena, and how does 
that differ from the functioning of a phase oscillator? 
Among others, the change of frequency (phase evolu- 
tion speed) should be emphasized. While any perturba- 
tion on a phase oscillator produces only a phase shift, 
the frequency in a high-dimensional oscillator can be 
modulated [l5|. When such oscillators are assembled 
to form a network, the resulting system can generate 
distributed frequencies states and lead to the temporal 
formation and breakdown of synchrony of the system. 
These behaviors are observed in many complex systems 
[1, 0, H, @, [13 , and in fact variable frequency effects have 
been noted in some them. Several theoretical studies 
have pointed to the variable frequency effect on complex 
systems [H [H [H [H . 

In the present paper, we introduce networks of in- 
teracting limit cycle oscillators, whose frequencies are 
modulated by changes in their amplitudes, and demon- 
strate that amplitudes of oscillators spontaneously take 
distributed values; i.e. the frequencies are distributed. 



Additionally, we observe that there are two types of be- 
haviors depending on network topologies: arrangement 
of frequencies at a regular interval and continuous irreg- 
ular modulation of frequencies. The analysis reveals the 
mechanism that generates distributed frequencies states. 
We also show how the topology of the networks regulates 
the behaviors of the system. 

Let us consider a limit cycle oscillator that is weakly 
stable in one direction transverse to the rotation direc- 
tion. If one produces perturbations of intensity compa- 
rable to the attraction of the limit cycle, the orbit widely 
deviates from the limit cycle, introducing considerable 
modulation into the frequency. In order to represent 
this effect in a simple model, we introduce a two- variable 
oscillator by assuming that the deviations in transverse 
directions are projected in a one-dimensional variable r 
called amplitude. Using as the phase variable defined 
along with the limit cycle, and taking the first order of r 
into account, we get the model equations. 



r = -T{e)r 



(perturbation). 



(1) 



Q.{9) represents the modulation of the frequency by the 
amplitude change, and T[9) indicates the decay rate of r 
which is assumed to be small. If the phase evolution is 
sufficiently faster than the change in the amplitude (weak 
stability assumption), the phase dependence of T[0) is 
normalized and replaced with the parameter 7. We adopt 
ijjQ = 1, Vt{0) = ijj ■ {1 + esin(27r6')) in this study. 

Next, we introduce the oscillator networks. Oscilla- 
tors are supposed to interact in a bidirectional way with 
an interaction matrix fey = kji G {0,1}. Note that the 
phase differences between interacting oscillators are not 
always small, and large differences in phases produce con- 
siderable perturbations in the amplitudes. We select the 
interaction function to agree with the diffusive coupling 
at the leading order, which satisfies dO / d/S.9\^r=Q 7^ 0, 
dr/dM\/s^r=o = 0, and d'^r / dM'^\i^r=o + 0, where M 
and Ar are phase differences and amplitude differences 



(C) 




^0.5 VrtiWvVvW^^^ 

■§0.4 ■ VWVJVVVVVIM^ / 

iO-2 

Co 

0.1 

^00 300 



400 time 500 



600 



700 




0.1 |AAAAAAAAAAAAAAAAAAAAAAAA/\AAAAAArtAAAAAAAAAAAAAAA/^rf 

^00 300 



400 time 500 



600 



700 



FIG. 1: (color online) Behaviors of coupled oscillator networks indicated by Fig. A and D are shown at Fig. B, C and E, F, 
respectively, e = 0.5, in B, C, E and F, and oj = 1 in C and F. B and E show the variances of amplitudes of oscillators in 
network A and D as functions of lo. Dashed lines indicate the slope of C and F show the time series of amplitudes of 

oscillators in network A and D. 



between the oscillators. Thus, the models of the oscilla- 
tor networks are given by 

Oa = I + uj {1 + esm{2n6a)) Ta 

+ DeJ2 ^0.0 sin {2Ti{ep ~ 9^)) , (2) 

Ta = -7rQ + A- {1 - cos(27r(6'/3 - 6*0))} ,(3) 
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where a, (3 are oscillator indexes. In this study, we set 
Dg = Dr = 0.01, 7 = 0.1, e = 0, 0.5 and changed the con- 
trol parameter tu = 0.2 ~ 2, which indicates the magni- 
tude of the modulation of the frequency by the amplitude 
change. All simulations run from uniformly distributed 
random initial conditions in Tq € (0, 1], 9a & (0, 1]. 

Now we take up two examples of networks indicated 
by Fig. [1]-A and Fig. [1]-D and demonstrate two distinct 
behaviors of the oscillator-networks. 

Figure [TJC shows the time evolution of amplitudes 
of oscillators in network A (Fig. [1]-A) after they reach 
the stationary states. The case with ui = 1, e = 0.5 
is shown. In the figure, amplitudes are distributed and 
remain around the discrete stationary points which are 
arranged at regular intervals of 0.1. The fluctuations of 
amplitudes around the stationary points are small and 
seem to be periodic or quasi-periodic. In other words, we 
clearly observe the layers of amplitudes. The variances of 
the amplitudes are plotted against to in Fig. [TJB. With 
uj < 0.38, all oscillators achieve complete synchroniza- 
tion (r^ = 0, 9a = 9p for all a,/3), and no fluctuation 
appears [l3|. With uj > 0.38, amphtudes are distributed 
and the averaged values of amplitudes are independent 
from UJ (data not shown). In the region, the variances of 
amplitudes decay with l/uj"^. 

On the other hand, network D (Fig. [1]-D) exhibits a 
different behavior as shown in Fig. [1]-E and -F. At the 



steady state, the layers of amplitudes are partially broken 
and the amplitudes of oscillators c, f, g and h show large 
fluctuations. The widths of their fluctuations arc greater 
than the range between layers, and the fluctuations de- 
cay slowly with the increase of uj. Other oscillators (a, 
b, d and e) still stay around the stationary points. We 
emphasize that network D is constructed by changing 
only one interaction from network A: c-d to c-g. Thus, 
behaviors are quite sensitive to the topology of networks. 

Summarizing the results, the amplitudes of oscillators 
in networks are widely distributed when the frequency 
modulation effect (uj) is large. We also found two distinct 
behaviors in amplitudes depending on the topologies of 
networks: (i) amplitudes of all oscillators in network form 
layers, and (ii) amplitudes of a part of oscillators show 
large fluctuations across the layers. Type (i) behavior is 
less frequently observed in larger networks. When 1000 
networks were generated with 10 and 50 oscillators, re- 
spectively, 468 networks with 10 oscillators showed type 
(i) behavior and only 25 networks with 50 oscillators 
showed type (i) behavior (iTj. Meanwhile, both widely 
fluctuated and less fluctuated oscillators always coexist 
in large networks. Figure [5] shows scatter plot distribu- 
tions of the averages of amplitudes and the variances of 
amplitudes for two random networks which have 100 and 
1000 oscillators. The variances are widely distributed. 
For less variant oscillators, the averaged amplitudes show 
the order in regular intervals as seen in Fig. [T]-C. 

In the following, we illustrate the mechanisms of gen- 
erating both behaviors. In the analysis e = is used for 
simplicity. 

We begin with investigating the mechanism of form- 
ing layers. In order for the amplitudes to take stationary 
values, interaction terms should be approximately con- 
stant. There are two possible mechanisms; one is phase 
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FIG. 2: (color online) Scatter plot distributions of averages 
and variances of amplitudes of all oscillators in two large net- 
works with 100 and 1000 oscillators. The networks are gen- 
erated randomly with connection probabilities 2/(number of 
oscillators), e — 0.5, uj = 1. Data of the 1000-osciIlator 
network are plotted with red open circles, and data of the 
100-oscillator network are plotted with blue filled squares. 



locking and the other is averaging. If the phase differ- 
ence between two interacting oscillators takes a constant 
value (a phase- locked state), the interaction term also 
becomes constant. However, it is unlikely that the phase 
differences between all interacting oscillators take con- 
stant values except in complete synchronization, where 
all oscillators assume an identical state (ra = 0,0a = Op 
for all a,f3). On the other hand, the averaging assump- 
tion for an interaction holds when the two oscillators have 
considerably different frequencies so that the time evolu- 
tion of amplitudes can not follow the alternation of the 
interaction term. Since the difference in amplitudes pro- 
duces a difference in frequencies, the averaging assump- 
tion is expected to hold. 

Assume that the speeds of the phase differences be- 
tween all coupled oscillators are sufficiently fast that in- 
teraction functions in cq. ([3]) arc approximately replaced 
with the averaged values (cos (27r(6'^ — Oa)) 0). Then, 
eq. ([3]) is solved and the approximate amplitude takes 
NaDj./j, where is the number of interactions for os- 
cillator a. Note that the approximate amplitudes are 
proportional to N^, and hence they take discrete values. 
With the parameters used in Fig. [T]-C, Dr/'y = 0.1 which 
agrees with the observed interval in the figure. 

Here we examine under which conditions the averag- 
ing assumption holds. The following equation expresses 
the phase difference between two interacting oscillators in 
which amplitudes are replaced with approximate values 



Map = [Na - Np) 



DrUJ 
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De sin {2TTA0ap) ■ (4) 



We also assume that the interaction effects from other 
oscillators are replaced by averaged values. U\Na—Np\ > 
and DrUj/{'yDg) ^ 1 hold, the interaction has little 
effect and a uniform oscillation with frequency F^p = 



Dr{Na — Np)u)/^ is an approximate solution of eq. ([4]). 
Putting the solution in cq Q, wc get the equation for 
the amplitude 

= -ITa + Dr{Na - 1) + A- {l " COS (27rA^„^) } , (5) 

where the effects from other interactions are again re- 
placed by averaged values. The solution is 



DrNa 



Dr 



7 ^{2nFap)^ + r 



:sm{2TTFapt + 0o), (6) 



where is the initial phase. The second term denotes 
that the variance of the amplitude decreases along with 
the increase of io/'y. This is also consistent with the de- 
crease of fluctuation in Fig. [T]-B, since 1 / V variance oc 
Fap oa Lu. Thus, the averaging assumption is fulfilled. 

When Na ~ Np holds, the averaging assumption docs 
not hold. The synchronous state {AOap = 0) becomes 
the attractor for eq. Thus, we get Ta — rp = 

(A'q, — l)Dr/7 from eq. ([5]). Although this means that the 
assumed relation for amplitudes (r^ =fp= NaDr/"f) in 
eq. (|4]) is violated, both amplitudes change simultane- 
ously and the relation = rp still holds. Hence, the 
synchrony of oscillators a and /3 is sustained as a special 
phase-locked case. 

To sum them up, two interacting oscillators with dif- 
ferent Na have sufficient difference in frequency to be 
averaged, and the interaction is approximately replaced 
by the averaged value. On the other hand, two oscillators 
with the same Na have the same frequency on average, 
and they synchronize with each other, an interaction that 
produces little effect on the amplitudes. Let us apply the 
analysis to network A and check the relation between the 
topology of the network and the observed layers in Fig. 
[T]-C. In network A, d-e and g-h pairs have the same num- 
bers of interactions {Nd = Ne = 2, Ng = Nh = 3); thus, 
they are supposed to be synchronizing pairs. Therefore, 
the averaged amplitudes of oscillators are given by 



7 



a 



a, ••■ 



(7) 



where N'a denotes the number of effective interactions 
after eliminating the interaction of synchronizing pairs. 
They show good agreement with the order of amplitudes 
in regular intervals as in Fig. [1]-C. Besides, the phase 
difference between oscillators d and e is always less than 
0.013 and oscillators g and h have exactly the same phase 
in the simulations. 

Now we investigate the continuous irregular motion of 
amplitudes observed in the network D. In short, it is 
caused by the indeterminacy of a consistent set of syn- 
chronizing pairs. This leads to continuous changes of the 
synchronizing pairs, and irregular motion arises. Here we 
illustrate these processes in detail using network D as an 
example. 



Following the above rules, oscillators / and g in net- 
work D synchronize {Nf = Ng = 4), and their ampli- 
tudes converge to ZDr/^- Additionally, oscillator c and 
/i, which interact with have the same frequency, since 
Nc = Nh = Ng — 1 = 3. Thus, g and either of two oscil- 
lators start to synchronize, which brings about another 
change in amplitudes of the synchronizing pair [c-g or g- 
h). Their amphtudes converge to 2Dr/^- However, this 
introduces a difference in amplitude between oscillator / 
and (7, i.e. a difference in frequencies, and hence their 
synchrony is disrupted. This effect also destroys the sec- 
ond synchrony {c-g or g-h). Thus, the system returns to 
the original no-synchrony state. In this way, the pairs of 
synchronization in the network D are never settled and 
the amplitudes of oscillators c, /, g and h continue to 
change spontaneously. For larger networks, it becomes 
more difficult to determine the pairs of synchronizing os- 
cillators in a consistent manner. Thus, complete layered 
behaviors becomes rare in larger networks. 

In summary, we have investigated networks of limit 
cycle oscillators the frequencies of which change widely 
along with the amplitude changes, and report novel com- 
plex behaviors: the layering and the continuous large 
fluctuation of amplitudes. In the analysis, we explained 
a mechanism that determines the behaviors of amplitudes 
from the local structure of the network. The behavior of 
amplitude primarily depends on the number of connec- 
tions the oscillator has, but also depends on the states of 
the connecting oscillators which are determined by their 
numbers of connections. Therefore, the global structure 
of the network topology affects the behavior of each os- 
cillator. 

The dependence of the frequency on amplitude changes 
is key for a variety of phenomena. In layered states, the 
differences in amplitudes originate from the differences 
in frequencies, and this assures the decoupling of oscil- 
lator interactions by averaging. Thus, the differences in 
amplitudes are sustained. The same effect enables the 
globally coupled oscillators to form stable multi-cluster 
states [ulllSl- By contrast, in fluctuating states, ampli- 
tude changes result in switching between decoupled and 
influential interactions. This emergence and annihilation 
of effective interactions sustainably drives the fluctuation 
of amplitudes. 

In this analysis, we also show that the topology of 
networks determines the behaviors of the oscillator net- 
works, and reveal how a small difference in topology 
brings about qualitative change in behaviors. This in- 
dicates that exploring the detailed topology of oscil- 
lator networks |14l | is as important as analyzing them 
statistically 0. 

This change in effective frequencies (time scales) and 
the associated change of effective interactions are impor- 
tant properties of complex systems. Such behaviors are 
frequently observed in biological systems as adaptations 



of time scales. For example, nerve systems are thought 
to code information by changing the firingfrequencics of 
neurons as well as their phase relations [7[. The com- 
plex dynamics of both properties are expected to work 
as unifled systems [1]. These relations are also observed 
in adaptive behaviors of single cell organisms which con- 
tains a variety of time scales 0, [l^l . Our study extracts 
the mechanism of changing time scales from a general 
limit cycle system, and reveals the origin of these com- 
plex phenomena. 

The author is grateful to K. Kaneko, S. Ishihara, and 
K. Fujimoto for helpful suggestions. 
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